Topological Aspect of Knotted Vortex Filaments in Excitable MediaEl 
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Scroll waves exist ubiquitously in three-dimensional excitable media. It's rotation center can 
be regarded as a topological object called vortex filament. In three-dimensional space, the vortex 
filaments usually form closed loops, and even linked and knotted. In this letter, we give a rigorous 
topological description of knotted vortex filaments. By using the (^-mapping topological current 
theory, we rewrite the topological current form of the charge density of vortex filaments and use this 
topological current we reveal that the Hopf invariant of vortex filaments is just the sum of the linking 
and self-linking numbers of the knotted vortex filaments. We think that the precise expression of 
the Hopf invariant may imply a new topological constraint on knotted vortex filaments. 
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Scroll waves exist ubiquitously in three-dimensional 
excitable media. They have been directly observed in 
three-dimensional chemical p]], physical @, Q, and biolog- 
ical systems Specific examples include the Belousov- 
Zhabotinskii reaction^, the cardiac muscleQ, and the 
oxidation of CO on platinum^. Recently, scroll waves 
have drawn great interest and have been studied inten- 
sively in many ways because it is believed to be the mech- 
anism of some re-entrant cardiac arrhythmias and fibril- 
lation which is the leading cause of death in the industri- 
alized world 0, H IE EE|. The investigation of properties 
of scroll waves in excitable media provides insights into 
the possible behaviors of these processes in the cardiac 
tissue of animals whose heart wall is thick enough for 
three-dimensional effects to be significant. 

On the other hand, the global behavior of a scroll wave 
is well described by the motion of its approximate ro- 
tation center which is a line defect known as a vortex 
filament and usually defined in terms of a phase sin- 
gularity. In three-dimensional excitable media, the vor- 
tex filaments is commonly a closed ring, and these vor- 
tex filaments can form linked and knotted rings which 
contract to compact, particle-like bundles [H EE EE 
EI EE EE El El EE HE EH- The existence of the 
knotted vortex filaments can be regarded as a topolog- 
ical phenomena in excitable media. In most previous 
works[IE El El EE EE EH, the topological arguments 
have been applied to help us understand the topologi- 
cal properties of the vortex filaments and some most im- 
portant topological constraints on behaviors of the vor- 
tex filaments have been investigated. These topologi- 
cal rules may have some important applications in prac- 
tice. In particular, the topological constraint on knot- 
ted vortex filaments is believed to relate to topological 
characteristic numbers of knotted vortex filament fam- 
ily, such as the winding, the self-linking and the linking 
numbers EE EE EH- So in research into knotted vortex 



filaments in excitable media, one should play much at- 
tention to these knot characteristics. In this letter, we 
will use the topological viewpoint to study the knotted 
vortex filaments with the Hopf invariant [23| which can be 
used to describe the linkages of knot family in mathemat- 
ics, and reveal the inner relationship between the Hopf 
invariant and the topological characteristic numbers of 
the knotted vortex filaments. We think that this inner 
relationship may imply a new topological constraint on 
knotted vortex filaments. 

We chose to work with a general two- variable reaction- 
diffusion system which mathematical description in terms 
of a nonlinear partial differential equation. This equation 
is written as 

dtu = f{u,v) + D u \7 2 u, d t v = g{u,v)+D v V 2 v, (1) 

where u and v represent the concentrations of the 
reagents; V 2 is the Laplacian operator in three- 
dimensional space; f(u,v) and g(u,v) are the reaction 
functions. Following the description in Ref.[l6[, we de- 
fine a complex function Z — (f) 1 + itfi 2 , where <f> 1 = u — u* 



and d> = v 



Here u* and v* are the concentrations 



of the vortex filaments. 

As pointed out in Ref . [24| , the sites of the vortex fil- 
aments are just the isolated zero lines of the complex 
function Z — (j) 1 +i(j> 2 in three-dimensional space. By us- 
ing the topological viewpoints, Zhang et al. Ell derived a 
topological expression of the charge density of these zero 
lines, i.e., the vortex filaments, which is written as 
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where Wi is the topological charge of the l-th vortex fil- 
ament L/.This expression reveals the topological struc- 
tures of vortex filaments. In our topological theory of 
knotted vortex filaments, the topological structure of the 
vortex filaments will play an essential role. In order to 
study the knotted vortex filaments more conveniently, we 
firstly rewrites the charge density of vortex filaments as 
a topological current. Now we begin to derive the topo- 
logical current form of the charge density of vortex fila- 
ments. We know that the complex function Z — (j) 1 + i(f> 2 
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can be regarded as the complex representation of a two- 
dimensional vector field Z = (cb 1 ,^ 2 ). Let us define the 



unit vector: n a = 



TOT 



(a =1,2; 



b a (f> a = Z*Z). It 



is easy to see that the zeros of Z are just the singulari- 
ties of n. Using this unit vector ft, we define a topological 
current 

j* = ^ k e ab d jn a d k n\ i,j,k= 1,2,3. (3) 
Apply the ^-mapping theory (22l. l25l, |26| , one can obtain 



where 



l f ijk 



lb ; 



(4) 



£abdj4> a dk(i> is the Jacobian vector. 



This delta function expression of the topological current 
j % tells us it doesn't vanish only when the vortex filaments 
exist, i.e., Z = 0. The sites of the vortex filaments deter- 
mine the nonzero solutions of j 1 . The implicit function 
theory shows that under the regular condition [27| 

^0, (5) 

x 

the general solutions of 

<f> 1 (x,t) = 0, cb 2 (x,t) = (6) 
can be expressed as 

x 1 =x}(t,s), x 2 = xj{t,s), x 3 =xf{t,s), l = l,2,-,N, 

which represent the world surface of N moving isolated 
vortex filaments Li(l — 1,2, ■ • ■ , N) with string parameter 
s. In delta function theory [28j, one can prove that in 
three-dimensional space, 



(8) 



where £>(£) = l^ k e mn ^^ and Ej is the l-th planar 
element transverse to L[ with local coordinates (w 1 ,?/ 2 ). 
The positive integer (3i is the Hopf index of ^-mapping, 
which means that when x covers the neighborhood of 
the zero point xi(s,t) once, the vector field </> covers the 
corresponding region in <fi space for /3j times. Meanwhile 
the direction vector of Li is given by[H, [2(| 



(9) 



dx* _D i (4>/x) 
H^ 1 ~ D{<p/u) 



Then considering Eqs.© and Eqs.©, we obtain the in- 
ner structure of j 1 , 



N 

^Pim / dx % 5 3 (x-xi), 

1 = 1 



(10) 



where rji = sgnD(f t ) = ±1 is the Brouwer degree of 
^-mapping, with r\i = 1 corresponding to the vortex fil- 
ament and r/i = — 1 corresponding to the antivortex fil- 
ament. Compare Eqs.([2]) and Eqs. fTU)) we see that the 
topological current j is just the charge density vector p of 
the vortex filament in Ref.[lB|. In our theory, the topo- 
logical charge of the vortex filament Li is 



= / j ■ da = Wi = 



(11) 



in which Wi is just the winding number of Z around 
Li , the above expression reveals distinctly that the topo- 
logical charge of vortex filament is not only the wind- 
ing number, but also expressed by the Hopf indices and 
Brouwer degrees. The topological inner structure showed 
in Eq. tfTTj) is more essential than in Eq.©, this is just the 
advantage of our topological description of the vortex fil- 
aments. 

Now let us begin to discuss the topological properties 
of knotted vortex filaments in excitable media. It is well 
know that the Hopf invariant is an important topological 
invariant to describe the topological characteristics of the 
knot family. In our topological theory of knotted vortex 
filaments, the Hopf invariant relates to the topological 
characteristics numbers of the knotted vortex filaments 
family. In a closed three-manifold M the Hopf invariant 
is defined asEJEl 
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Aifd 6 x, 



(12) 



in which is a "induced Abelian gauge potential" con- 
structed with the complex function Z. The relation- 
ship between the Abelian gauge field strength Fij — 
diAj — djAi and the topological current j 1 is [261] 



^ ijk Fjk = ^ k e ab d jn a d k 
Air 2tt 



.1 : — <" = — <'■'"<>, i.<>,ir<>i,n . (13) 
Substituting Eq. (fTU|) into Eq. (fl"2")l . one can obtain 
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A,dx\ 



(14) 



It can be seen that when these N vortex filaments are N 
closed curves, i.e., a family of N knots £/(/ = 1, 2, N), 
Eq.CSD leads to 



H 



i N r 

i=i j & 



Ajdx % 



(15) 



This is a very important expression. Consider a trans- 
formation of complex function Z = e 1,6 Z, this gives the 
U(l) gauge transformation of Ai : A i = Ai + di6, where 
9 6 R is a phase factor denoting the U(l) gauge trans- 
formation. It is seen that the diQ term in Eq. (fl5[) con- 
tributes nothing to the integral H when the vortex fila- 
ments are closed, hence the expression (|15[) is invariant 
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under the U(l) gauge transformation. As pointed out in 
Ref. [TH , a singular vortex filament is either closed ring 
or infinite curve, therefore we conclude that the Hopf 
invariant is a spontaneous topological invariant for the 
vortex filaments in excitable media. 

It is well known that many important topological num- 
bers are related to a knot family such as the self-linking 
number and Gauss linking number. In order to discuss 
these topological numbers of knotted vortex filaments, 
we define Gauss mapping: 



: S 1 x S 1 -> S 2 , 



where m is a unit vector 



m(x,y) 



y-x 
\y- x\ 



(16) 



(17) 



where x and y are two points, respectively, on the knots 
and (in particular, when x and y are the same point 
on the same knot £ , ft is just the unit tangent vector T 
of £ at x ). Therefore, when x and y , respectively, cover 
the closed curves £fc and £; once, n becomes the section 
of sphere bundle S 2 . So, on this S 2 we can define the 
two-dimensional unit vector e = e(x, y). e, rh are normal 
to each other, i.e. , 

e l • & 2 = e i • rh = e*2 • m = 0, 

ei • ei = e*2 • e*2 = m • rh = 1. (18) 

In fact, the gauge field A can be decomposed in terms of 
this two-dimensional unit vector e: Ai = e a be a die b — di9, 
where is a phase factor (2f| Hfjj]. Because the (di9) term 
does not contribute to the integral H, Ai can in fact be 
expressed as 

Ai = e ab efd i e h . (19) 
Substituting it into Eq.(14), one can obtain 



1 N f 
H = — 22 W k f £abe a (x, y)die b (x, yjda 

n k=i ^ k 



(20) 



Noticing the symmetry between the points x and y in 
Eq. (fT7|) . Eq. ([20|) should be reexpressed as 



N 

H =^~Y W kWi l l eabd^d^dx 1 Ady j . (21) 



fcj=i 



In this expression there are three cases: (1) £/. and £z 
are two different vortex filaments (£& ^ and a? and y 
are therefore two different points (x ^ y); (2) and 
are the same vortex filaments (£& = but x and y are 
two different points (x ^ y)\ (3) and £j are the same 
vortex filaments = and x and ?/" are the same 
points (x = y). Thus, Eq. (|2~lj) can be written as three 



terms: 



H 



N I 

k=l(k=l,x"=jiy) 
N 



iff e ab d i e a d j e b dx i A dy j 



Y^W 2 ! e ab e a d ie b dx l 

1 1 J £ h 



fc=l 

N 



y 7r WkWi f f e - bd ^ aa 



e b dx t A dy J . (22) 



By making use of the relation e a i,die a dje b = hfh ■ (dim x 
djrh), the Eq.(|22J) is just 



2 f d> m*(dS) 



1 * / 

+ YY, W H ^ b e a d ie b dx 



fc=i 



+ y ^ WkWi <f i ™*( ds ), ( 23 ) 

where rh*(dS) = rh- (difh x djrnjdx" 1 Ady J (x ^ y) denotes 
the pullback of the S 2 surface element. 

In the following we will investigate the three terms in 
the Eq.([23]l in detail. Firstly, the first term of Eq.([23]) is 
just related to the writhing number Wr(^) of ft.|29l [30|| 



47T 



rh*(dS). 



(24) 



For the second term, one can prove that it is related to 
the twisting number Tw(^k) of 



1 

2^ 



e ab e a d z e b dx l 



1 

27 



(T x V) ■ ^ 



(25) 



where T is the unit tangent vector of knot at x (m = T 
when x = y) and K is defined as e a = e ah l^ h (T7 iT,e = 
f X F). In terms of the White formula^ ^ 



(26) 



Si(^) = Wr(&) + Ti«(^), 



we see that the first and the second terms of Eq. (|23]) just 
compose the self-linking numbers of knots. 

Secondly, for the third term, one can prove that 



1 



-in 



6= J£ 



— d> 6 m*(dS) 



= — e ^ k i dx l I dv j ^ ~ ^ 



(27) 
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where is the Gauss linking number between 

and & [3l]]. Therefore, from Eqs.(El, J35J), J2BJ) and ([2T|. 
we obtain the important result: 

N N 

H = Y J WlSL{^ k )+ Y, WkWtLk^b). (28) 
fc=i fc,;=i(fc^i) 

This precise expression just reveals the relationship be- 
tween H and the self-linking and the linking numbers of 
the vortex filaments knots family. Since the self-linking 
and the linking numbers are both the invariant char- 
acteristic numbers of the vortex filaments knots family 
in topology, H is an important topological invariant re- 
quired to describe the linked vortex filaments in excitable 
media. 

So far, we obtained a more essential topological for- 
mulary of charge density of knotted vortex filaments and 
revealed the topological inner relationship of between the 
Hopf invariant and the self-linking and the linking num- 
bers of knotted vortex filaments. We conclude that the 



Hopf invariant is just the sum of the linking and self- 
linking numbers of the knotted vortex filaments. In the 
present study, it should be pointed out that, when we dis- 
cussed the topological properties of the vortex filaments, 
the regular condition D l (&) ^ must be satisfied. Now 
the question is coming, when this condition fails, what 
will happen about the vortex filaments? The answer is 
related to the evolution of the vortex filaments [2a, |26| |. 
As we all known that the evolution of the vortex fila- 
ments which include generating, annihilating, splitting, 
or merging may obey some topological constraints. These 
constraints usually relate to the topological numbers of 
the vortex filaments, such as the topological charges and 
linking numbers. So it is naturally to think that the pre- 
cise expression Eq. (|28p of Hopf invariant may imply a 
new topological constraint on the behavior of the knot- 
ted vortex filaments. What the rigorous description of 
this possible new constraint is will be investigated in our 
further works. 
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